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A  most  intuitive  realization  of  a  qubit  is  a  sin¬ 
gle  electron  charge  sitting  at  two  well-defined  po¬ 
sitions,  such  as  the  left  and  right  sides  of  a  double 
quantum  dot.  This  qubit  is  not  just  simple  but 
also  has  the  potential  for  high-speed  operation, 
because  of  the  strong  coupling  of  electric  fields  to 
the  electron.  However,  charge  noise  also  couples 
strongly  to  this  qubit,  resulting  in  rapid  dephas¬ 
ing  at  nearly  all  operating  points,  with  the  ex¬ 
ception  of  one  special  “sweet  spot.”  Fast  dc  volt¬ 
age  pulses  have  been  used  to  manipulate  semi¬ 
conductor  charge  qubits  mm,  but  these  previ¬ 
ous  experiments  did  not  achieve  high-fidelity  con¬ 
trol,  because  dc  gating  requires  excursions  away 
from  the  sweet  spot.  Here,  by  using  resonant  ac 
microwave  driving,  we  achieve  coherent  manipu¬ 
lation  of  a  semiconductor  charge  qubit,  demon¬ 
strating  an  Rabi  frequency  of  up  to  2  GHz,  a 
value  approaching  the  intrinsic  qubit  frequency 
of  4.5  GHz.  Z-axis  rotations  of  the  qubit  are 
well-protected  at  the  sweet  spot,  and  by  using  ac 
gating,  we  demonstrate  the  same  protection  for 
rotations  about  arbitrary  axes  in  the  X-Y  plane 
of  the  qubit  Bloch  sphere.  We  characterize  oper¬ 
ations  on  the  qubit  using  two  independent  tomo¬ 
graphic  approaches:  standard  process  tomogra¬ 
phy  [8|  [9j  and  a  newly  developed  method  known 
as  gate  set  tomography  m •  Both  approaches 
show  that  this  qubit  can  be  operated  with  pro¬ 
cess  fidelities  greater  than  86%  with  respect  to  a 
universal  set  of  unitary  single-qubit  operations. 

Coherent  control  of  qubits  with  resonant  microwaves 
plays  an  essential  role  in  realizing  precise  single  mmm 
and  two  m  qubit  gates  in  solid  state  quantum  com¬ 
puting  architectures  m-  In  electrically-controlled  quan¬ 
tum  dots,  driven  coherent  oscillations  have  been  demon¬ 
strated  in  spin-based  qubits  using  electron  spin  reso¬ 
nance  pn,  electric-dipole  spin  resonance  [2QH23] .  or 
resonant  exchange  gates  [24],  with  typical  rotation  rates 
(Rabi  frequencies)  on  the  order  of  1  to  100  MHz.  Here  we 
demonstrate  fast  coherent  operation  and  full  control  on 
the  Bloch  sphere  of  a  charge  qubit  in  a  double  quantum 
dot  in  a  silicon/ silicon- germanium  (Si/SiGe)  heterostruc- 
ture. 

The  charge  qubit  we  study  is  formed  by  three  electrons 
in  a  Si/SiGe  double  quantum  dot  (Fig.  la)  [7,  11].  Fig.  lc 
shows  the  qubit  energy  level  diagram  as  a  function  of 
the  detuning  6,  which  is  controlled  by  gate  GL  or  GR 


(Fig.  la).  The  states  |2, 1)  =  | L)  and  1 1,  2)  =  \R)  are  the 
ground  states  of  the  system  at  negative  and  positive  6, 
respectively,  and  these  states  anticross  at  e  ==  0.  A  state 
|R)e,  of  higher  energy  (>14  GHz  with  respect  to  1 1)  near 
6  =  0  [12])  and  not  visible  in  the  energy  range  shown  in 
the  figure,  affects  the  dispersion  of  energy  levels  shown, 
but  is  otherwise  unimportant  because  its  occupation  is 
negligible  for  the  resonant  driving  demonstrated  here. 

We  focus  on  detunings  near  6  ~  0,  where  an  avoided 
crossing  is  formed  between  states  | L)  and  | R)  with  tun¬ 
nel  coupling  strength  Ai.  The  logical  qubit  states  are  the 
energy  eigenstates  at  6  =  0:  |0)  =  (\L)  +  \R))/y/2  and 
|1)  =  (| L)  —  \R))/y/2.  These  states  can  be  adiabatically 
evolved  to  states  | L)  and  \R)  at  moderate  negative  6, 
enabling  robust  readout  of  Pi ,  the  probability  of  occupa¬ 
tion  of  |1)  at  the  end  of  qubit  manipulation,  performed  by 
monitoring  the  current  /qpc  through  the  charge-sensing 
quantum  point  contact  (Fig.  la). 

We  first  characterize  the  system  spectroscopically.  As 
shown  schematically  by  the  black  solid  line  in  Fig.  lc, 
state  | L)  is  first  prepared  by  waiting  longer  than  the 
charge  relaxation  time  (7\  =23.5  ns)  at  a  detuning  ap¬ 
propriate  for  initialization  and  readout,  er  ~  —160  fieV 
lT2l .  We  then  ramp  the  gate  voltage  over  a  time  of  4  ns  to 
change  the  detuning  to  a  value  near  6  =  0,  adiabatically 
evolving  state  | L)  to  state  |0),  completing  the  qubit  ini¬ 
tialization.  At  this  detuning,  a  10  ns  microwave  burst  is 
then  applied  to  gate  GR.  When  the  microwave  frequency 
is  resonant  with  the  splitting  between  the  qubit  energy 
levels,  excitation  occurs  from  |0)  to  |1).  The  resulting 
probability  of  state  |1)  is  measured  by  ramping  the  de¬ 
tuning  adiabatically  over  approximately  2  ns  back  to  6r, 
transforming  |0)  to  |  L)  and  |1)  to  |  R),  and  measuring 
the  change  in  /qpc  Q2].  The  details  of  the  measure¬ 
ment  procedure  and  conversion  to  probability  are  pre¬ 
sented  in  the  Supplementary  Information  SI.  Fig.  lb 
shows  the  resulting  spectroscopy  of  the  qubit  energy  lev¬ 
els.  We  find  a  good  agreement  between  the  spectroscopic 
measurement  with  the  calculated  lowest-energy  excita¬ 
tion  spectrum  (green  dashed  curve)  with  Hamiltonian 
parameters  measured  in  Ref.  m-  Because  the  low-lying 
excited  state  on  the  right  dot  |  R)e  affects  energy  disper¬ 
sion,  we  measure  a  minimum  resonant  frequency  of  ap¬ 
proximately  4.5  GHz  near  6  =  0,  typically  lower  than  the 
tunnel  coupling  2Ai /h  =  5.2  GHz  used  for  calculation. 
At  this  sweet  spot  dE/ds  =  0,  and  the  energy  levels  are 
first  order  insensitive  to  detuning  noise  00  El]. 

Coherent  oscillations  between  qubit  states  |0)  and  |1) 
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Figure  1.  Si/SiGe  quantum  dot  device,  qubit  spectroscopy,  and  coherent  Rabi  oscillation  measurements,  a, 

SEM  image  of  a  device  lithographically  identical  to  the  one  used  in  the  experiment,  with  the  locations  of  the  double  dot 
indicated  by  white  dashed  circles  HU.  The  current  through  the  quantum  point  contact  (QPC)  /qpc  is  used  for  charge  sensing 
via  a  measurement  of  its  change  when  microwave  component  is  added  to  a  voltage  sequence  applied  to  gate  GR.  b-c,  Qubit 
energy  levels  and  microwave  spectroscopy,  b,  Probability  Pi  of  the  state  to  be  |1)  at  the  end  of  the  driving  sequence  as  a 
function  of  detuning  s  and  excitation  frequency  /ex  of  the  microwave  applied  to  gate  GR.  For  this  experiment,  the  repetition 
rate  of  the  driving  sequence  is  ~  15  MHz  and  the  microwave  pulses  are  of  duration  10  ns.  The  base  value  of  the  detuning 
sr  used  for  both  readout  and  initialization  is  e  =  —160  fieV .  Pi  is  large  when  the  microwave  pulse  is  resonant,  so  that  it 
excites  the  system  from  the  ground  state  to  the  excited  state.  Dashed  green  curve  shows  fit  to  calculated  energy  difference 
between  ground  state  and  lowest-energy  excited  state  of  the  three-level  model  of  Ref.  [12]  (see  also  Supplementary  Information 
S2).  The  third  level  of  higher  energy  affects  the  dispersion  of  energy  levels  shown,  but  is  otherwise  unimportant  because  its 
occupation  is  negligible  for  all  experiments  shown,  c,  Diagram  of  the  calculated  energy  levels  E  versus  detuning  e,  including 
the  ground  states  of  the  (2, 1)  and  (1,2)  charge  configuration,  | L)  and  | R)  respectively,  and  logical  states  |0)  =  (| L)  +  \R))/y/2 
and  |1)  =  (| L)  —  \R))/y/2.  Black  solid  line  inset:  pulse  sequence  used  for  Rabi  oscillation  and  spectroscopy  measurements,  d-e, 
Rabi  oscillation,  d,  Pi  as  a  function  of  the  voltage  Vgl  and  microwave  pulse  duration  tb  with  /ex  =  4.54  GHz  and  excitation 
amplitude  Vac  =  70  mV.  e,  Line-cut  of  Pi  near  Vgl  =-390  mV,  showing  «  1  GHz  coherent  Rabi  oscillations.  Red  solid  curve 
shows  a  fit  to  an  exponentially  damped  sine  wave  with  best  fit  parameter  TJ  =  1.5  ns.  f-g,  Dependence  of  the  Rabi  oscillation 
frequency  on  the  microwave  amplitude,  f,  Pi  as  a  function  of  Vac  and  tb  with  /ex  =  4.54  GHz,  which  demonstrates  that  the 
Rabi  oscillation  frequency  /Rabi  increases  as  the  amplitude  of  the  microwave  driving  is  increased,  g,  Rabi  oscillation  frequency 
/Rabi  as  a  function  of  Vac  with  fixed  /ex  =  4.54  GHz.  The  good  agreement  of  a  linear  fit  (red  dashed  line)  to  the  data  is 
strong  evidence  that  the  measured  oscillations  are  indeed  Rabi  oscillations,  with  the  Rabi  frequency  proportional  to  the  driving 
amplitude. 


are  implemented  by  applying  the  microwave  sequence 
shown  in  Fig.  lc,  using  an  ac  excitation  frequency 
/ex  ^4.5  GHz,  which  is  resonant  with  the  qubit  at  the 
sweet  spot,  and  an  ac  microwave  amplitude  Vac  =  70  mV 
(as  measured  at  at  the  arbitrary  waveform  generator;  the 
amplitude  at  the  sample  is  lower  because  of  attenuation 
of  the  coaxial  cables  in  the  dilution  refrigerator).  Fig.  Id 
shows  the  resulting  microwave- driven  Rabi  oscillations  in 
Pi ,  measured  as  a  function  of  the  microwave  burst  dura¬ 
tion  tb  and  the  gate  voltage  Vgl,  which  determines  the 


base  level  of  6.  Fig.  le  shows  a  line  cut  at  Vl  =  —390  mV, 
corresponding  to  6  «  0,  revealing  periodic  oscillations  in 
Pi  at  a  frequency  /Rabi  ~  1  GHz.  The  dependence  of  Pi 
on  tb  is  well  fit  by  an  exponentially  damped  sine  wave 
(Fig.  le,  red  solid  curve),  with  the  best  fit  yielding  a 
coherence  time  X2*  =  1.5  ns. 

Fig.  If  reports  the  power  dependence  of  the  qubit  os¬ 
cillations,  showing  Pi  as  a  function  of  tb  and  the  peak-to- 
peak  microwave  amplitude  Vac.  The  observed  oscillation 
frequency  varies  linearly  with  Vaci  as  expected  for  Rabi 
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Figure  2.  Ramsey  fringes  and  demonstration  of  three  axes  control  of  ac-gated  charge  qubit,  a,  Schematic  of  the 
pulse  sequences  used  to  perform  universal  control  of  the  qubit.  Both  the  delay  te  and  the  phase  0  of  the  second  microwave 
pulse  are  varied  in  the  experiment,  b-c,  Experimental  measurement  of  Z-axis  rotation.  In  a  Ramsey  fringe  (Z-axis  rotation) 
measurement,  the  first  Xn/2  gate  rotates  the  Bloch  vector  onto  the  X-Y  plane,  and  the  second  Xn/2  gate  (0  =  0)  is  delayed 
with  respect  to  the  first  gate  by  te,  during  which  time  the  state  evolves  freely  around  the  Z-axis  of  the  Bloch  sphere,  b,  Pi  as  a 
function  of  Vgl  and  te  for  states  initialized  near  |Y).  c,  Line-cut  of  Pi  near  Vgl  =  -390  mV,  showing  «  4.5  GHz  Ramsey  fringe. 
Red  solid  curve  shows  a  fit  to  exponentially  damped  sine  wave  with  best  fit  parameter  T2  =1.3  ns.  d-e,  Two  axis  control  of 
the  qubit  d,  Pi  as  a  function  of  4>  and  te.  e,  Line-cut  of  d  at  0  =  0  (X-axis,  black),  90° (Y-axis,  red),  and  180° (-X-axis,  blue). 
The  coherent  z-axis  rotation  along  with  the  rotation  axis  control  with  0  demonstrates  full  control  of  the  qubit  states  around 
three  orthogonal  axes  on  the  Bloch  sphere. 


oscillations.  Taking  the  filtering  and  attenuation  in  the 
dilution  refrigerator  into  account,  we  estimate  the  peak- 
to-peak  amplitude  at  the  sample  for  Vac  =  120  mV  to 
be  1.1  mV.  The  maximum  measured  Rabi  frequency  of 
2  GHz  approaches  the  intrinsic  frequency  of  the  qubit 
(4.5  GHz)  and  is  observed  with  only  moderate  applied 
power.  These  oscillations  correspond  to  X-rotations  of 
the  qubit  on  the  Bloch  sphere,  and  for  the  fastest  rota¬ 
tions  observed  an  X^/2  gate  requires  a  time  of  125  ps. 

We  now  demonstrate  Z-axis  rotations  of  the  qubit 
through  the  performance  of  a  Ramsey  fringe  experiment, 
using  the  microwave  pulse  sequence  shown  schematically 
in  Fig.  2a.  We  first  prepare  the  state  | Y)  =  ^/l/2(|0)  + 
i|l))  by  initializing  to  |L),  adiabatically  changing  the  de¬ 
tuning  to  6  =  0  fieV  to  evolve  the  state  to  |0),  and  then 
performing  an  X^/2  rotation  at  5  =  0  /xeV.  Z-axis  rota¬ 
tion  results  from  the  evolution  of  a  relative  phase  between 
states  |0)  and  1 1) ,  given  by  (p  =  —te2Ai/h,  where  te  is 
the  time  spent  at  5  =  0.  The  resulting  state  is  rotated 
by  a  second  X^/2  microwave  pulse,  and  the  final  prob¬ 
ability  Pi  again  is  measured  by  adiabatically  projecting 


state  |1)  to  |P)  at  the  readout  position  eT.  Figs.  2b  and  c 
show  the  resulting  quantum  oscillations  of  the  qubit  state 
around  the  Z-axis  of  the  Bloch  sphere  and  correspond  to 
a  Ramsey  fringe  experiment.  By  fitting  the  oscillations 
to  an  exponentially  damped  sine  wave  (red  solid  curve), 
we  extract  a  dephasing  time  T2*  =  1.3  ns  and  an  oscil¬ 
lation  frequency  2 A /h  =  4.5  GHz,  the  latter  of  which  is 
consistent  with  the  spectroscopic  measurements  shown  in 
Fig.  lb. 

A  key  advantage  of  ac  gating  over  dc  pulsed-gating  is 
the  ability  to  choose  the  rotation  axis  to  point  in  an  ar¬ 
bitrary  direction  in  the  X-Y  plane  of  the  Bloch  sphere  by 
controlling  the  phase  of  the  applied  microwave  burst.  On 
resonance  in  the  rotating  frame,  the  Hamilton  takes  the 
form  H  =  cos ( 0)crx+ sin (0)cry,  where  the  cq  are  Pauli  ma¬ 
trices  and  (p  is  the  relative  phase  of  the  microwave  burst 
with  respect  to  the  first  pulse  incident  on  the  qubit  |25] . 
Controlling  the  phase  of  pulses  subsequent  to  the  initial 
pulse,  as  shown  in  Fig.  2a,  thus  allows  rotations  of  the 
qubit  around  any  axis  lying  in  the  X-Y  plane  of  the  Bloch 
sphere.  Fig.  2d  shows  measurement  of  Pi  using  such  a 
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Figure  3.  Hahn  echo  measurement,  a,  Schematic  pulse 
sequence  for  the  measurement  of  Hahn  echo  that  corrects 
for  noise  that  is  static  on  the  time  scale  of  the  pulse  se¬ 
quence  [M  H>,  27  .  b-c,  Typical  echo  measurement  with  fixed 
total  evolution  time  r  =  1.4  ns.  b,  Pi  as  a  function  of  Vgl 
and  delay  time  St  of  the  pulse  with  decoupling  Xn  gate. 
The  effects  of  static  inhomogeneities  are  minimized  at  St  =  0, 
and  oscillations  of  Pi  as  function  of  St  at  twice  the  Ram¬ 
sey  frequency  decay  with  St  at  the  inhomogeneous  decay  rate 
1/T2*.  The  magnitude  of  the  signal  at  St  =  0  as  the  wait  time 
r  is  varied  decays  at  the  homogeneous  decay  rate  I/T2.  c, 
Line  cut  of  Pi  near  Vgl  =  -390  mV  showing  oscillatory  signal 
at  twice  the  Ramsey  frequency,  ~  9  GHz.  Solid  red  curve  is 
a  fit  to  gaussian  envelope  with  fixed  T2*  =  1.3  ns,  consistent 
with  the  Ramsey  fringe  measurement,  d,  Echo  amplitude  as 
a  function  of  r.  The  solid  red  curve  is  a  gaussian  fit  with 
T^echo  =  2.1  ns.  Applying  the  Hahn  echo  sequence  increases 
the  dephasing  time,  indicating  that  a  significant  component 
of  the  dephasing  arises  from  low-frequency  noise  processes. 


pulse  sequence  as  a  function  of  both  0  and  te,  demon¬ 
strating  smooth  variation  in  Pi  to  changes  in  the  control 
parameters.  Fig.  2e  shows  line  cuts  of  Pi  at  0=0,  90°, 
and  180°  (corresponding  to  the  second  pulse  inducing  a 
7t/2  rotation  around  the  X,  V,  and  —X  axes,  respec¬ 
tively);  microwave  phase  control  clearly  enables  control 
of  the  phase  of  the  resulting  Ramsey  fringes.  Taken  to¬ 
gether,  the  data  summarized  in  Figs.  1  and  2  demon¬ 
strate  control  of  the  qubit  over  the  entire  Bloch  sphere, 
including  both  Z-axis  rotations  and  rotations  about  any 
arbitrary  axis  lying  in  the  X-Y  plane. 

We  characterized  decoherence  times  by  implementing 
a  Hahn  echo  mi  na  mi  of  the  ac-gated  charge  qubit  by 
applying  the  pulse  sequence  shown  in  Fig.  [3^i.  Inserting 
an  Xn  pulse  between  state  initialization  and  measure¬ 
ment,  which  is  performed  with  X^/2  gates,  corrects  for 
noise  that  is  static  on  the  time  scale  of  the  pulse  se¬ 
quence.  In  Fig.  §>  and  c,  while  keeping  the  total  free 
evolution  time  r  fixed,  we  sweep  the  position  of  the  de¬ 


coupling  Xn  pulse  to  reveal  an  echo  envelope  [3127].  The 
maximum  amplitude  of  the  observed  envelope  reveals  the 
extent  to  which  the  state  has  dephased  during  the  free 
evolution  time  r,  characterized  by  T2,  whereas  the  ampli¬ 
tude  decays  as  a  function  of  St  with  inhomogeneous  decay 
time  T2  .  The  oscillations  of  Pi  in  Fig.  If  and  c  are  ob¬ 
served  as  a  function  of  St  at  twice  the  Ramsey  frequency 
(2/Ramsey  ~  9  GHz)  and  are  well-fit  by  a  Gaussian  decay 
(red  solid  curve).  Fig.  shows  the  echo  amplitude  decay 
as  a  function  of  r,  where  for  each  r  the  echo  amplitude  is 
determined  by  fitting  the  echo  envelop  to  Gaussian  decay 
similar  to  Fig.  If-  By  fitting  the  echo  amplitude  decays 
to  a  Gaussian,  we  obtain  the  dephasing  time  T2  «  2.2  ns. 

The  ultimate  test  of  experimental  qubit  control  is  the 
demonstration  of  repeatable  quantum  logic  gates.  Al¬ 
though  7r/2  rotations  that  generate  the  Clifford  group 
are  commonly  demonstrated,  ac  control  allows  direct 
implementation  of  any  unitary.  We  therefore  demon¬ 
strated  (and  validated  with  quantum  tomography)  two 
distinct  gatesets:  (1)  high-fidelity  approximations  to 
{X^Z^};  and  (2)  a  set  of  three  arbitrarily  chosen 
near- unit  ary  operations  {Gi,  G2,  G3}.  We  used  standard 
quantum  process  tomography  (QPT)  [9,  [28]  to  charac¬ 
terize  the  first  gateset.  Fig.  [4](a)  shows  the  resulting  pro¬ 
cess  matrices  (x)  expressed  in  the  Pauli  basis:  solid  bars 
represent  the  ideal  “target”  quantum  processes,  while 
open  circles  show  the  results  of  QPT.  The  process  fi¬ 
delities  29  between  the  QPT  estimates  and  the  targets 
are  F  =  86%  and  F  =  91%  for  the  X^/2  and  Zn/2  opera¬ 
tions,  respectively.  QPT  analysis  is  not  entirely  reliable, 
because  QPT  relies  on  prior  knowledge  of  input  states 
and  final  measurements  that  are  implemented  using  the 
same  logic  gates  that  we  seek  to  characterize  (see  Sup¬ 
plementary  Information  S3).  Therefore,  we  also  applied 
a  technique  called  gate  set  tomography  (GST)  |T0j  that 
avoids  these  assumptions.  GST  characterizes  logic  gates 
and  state  preparation/measurement  simultaneously  and 
self-consistently,  by  representing  all  of  them  as  unknown 
process  matrices.  Because  this  frees  us  from  any  obli¬ 
gation  to  use  carefully  calibrated  operations,  we  applied 
GST  to  a  set  of  three  repeatable  but  uncalibrated  logic 
gates  that  we  denote  {Gi,G2,Gs}. 

Data  for  GST  are  obtained  from  many  repetitions  of 
several  specific  experiments,  each  described  by  a  specific 
sequence  of  operations:  (i)  initialize  the  qubit  in  state 
p;  (ii)  perform  a  sequence  of  L  £  [0 ...  32]  operations 
chosen  from  {Gi,  G2,  G3};  (ii)  perform  measurement  M. 
Statistical  analysis  (a  variant  of  maximum  likelihood  es¬ 
timation)  is  used  to  find  the  estimates  {G 1,  G2,  G3,  p,  A4} 
that  are  most  consistent  with  the  measurements  (see  the 
Supplementary  Information  S3  for  more  details).  Since 
we  did  not  set  out  to  implement  any  particular  rotations, 
we  compute  ex  post  facto  the  closest  unitary  rotations  to 
these  estimates,  and  define  these  as  the  “target”  gates. 
The  results  are  shown  in  Figs|4]>c.  Figj4]3  shows  the  el¬ 
ements  of  the  process  (x)  matrix  for  the  GST  estimates 
(triangles)  and  the  closest-unitary  “targets”  (solid  bars). 
Figj4f:  portrays  those  same  unitaries  as  rotations  of  the 
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Figure  4.  Quantum  process  tomography  and  gate  set  tomography  of  the  ac-gated  charge  qubit,  a,  Real  and 
imaginary  parts  of  the  elements  of  the  process  matrix  y  [28]  in  the  Pauli  basis  {I,X,  Y,  Z}  for  Xn/2  and  Zn/ 2  processes  : 
ideal  “targets”  (solid  bars),  and  standard  QPT  estimates  (open  circles),  b,  y  for  the  uncalibrated  operations  Gi,  G2,  and 
G3,  obtained  by  gate  set  tomography  1  Oj  (GST,  triangles)  and  standard  quantum  process  tomography  [9]  1121  [28]  (QPT,  open 
circles),  compared  to  target  gates  Ti,  T2,  and  T3  (solid  bars).  Since  these  gates  are  not  precalibrated,  the  target  gates  are 
defined  to  be  the  unitary  processes  closest  to  the  GST  estimates  of  G 1,  G2,  and  G3  in  Frobenius  norm.  GST  self-consistently 
determines  the  state  preparation,  gate  operations,  and  measurement  processes  m-  c,  Rotation  axes  on  the  Bloch  sphere, 
rotation  angle  0 ,  process  fidelities  obtained  by  GST  (XPjgst),  and  QPT  (XPjqpt)  for  three  processes  G 1,  G2,  and  G3.  Here,  the 
rotation  axis  and  angle  correspond  to  the  closest  unitary  operations  to  the  GST  estimate  (Ti,  T2,  and  T3);  the  process  fidelities 
are  also  taken  between  the  estimates  and  these  target  processes.  The  error  on  XP)qpt  was  estimated  by  repeating  QPT  using 
10  distinct  sets  of  input  and  output  states;  standard  deviations  are  reported.  GST  and  QPT  yield  consistent  results,  with 
process  fidelities  >  86%  for  all  gates. 
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Bloch  sphere. 

We  also  performed  QPT  on  the  G\  and  G3  gates  to 
confirm  and  validate  the  GST  results.  QPT  depends 
critically  on  known  gates  (used  to  prepare  known  input 
states  and  perform  known  measurements),  so  we  used 
the  GST  closest-unitary  approximations  for  G\  and  G3 
to  model  the  preparation  of  input  states  and  final  mea¬ 
surements  for  QPT  (G2  was  neither  necessary  nor  used). 
The  results,  shown  as  open  circles  in  in  Figj4]3,  are  con¬ 
sistent  with  the  GST  results.  To  distill  a  single  figure 
of  merit,  we  computed  the  process  fidelity  F  between 
our  estimates  (both  GST  and  QPT)  and  the  closest  uni¬ 
tary  rotation.  Its  interpretation  is  slightly  different  in 
this  context;  since  the  “targets”  were  computed  ex  post 
facto ,  1  —  F  quantifies  the  amount  of  incoherent  error 
in  our  gate  implementation  (whereas  with  a  pre-existing 
target,  it  also  quantifies  coherent  under/over-rotation  er¬ 
rors).  These  process  fidelities,  shown  for  both  GST  and 
QPT  in  Figgjy  are  consistent  both  with  each  other  and 
the  process  fidelities  calculated  in  Figj4^i. 

Coherent  microwave  ac-gating  of  a  semiconductor 
quantum  dot  charge  qubit  offers  fast  (  >GHz)  manip¬ 
ulation  rates  for  all  elementary  rotation  axes.  Moreover, 
because  all  gates  can  be  performed  at  the  sweet  spot 


where  the  decoherence  time  is  in  the  order  of  nanoseconds 
instead  of  ^  100  ps  isiei,  rotations  around  three  orthogo¬ 
nal  axes  of  the  Bloch  sphere  with  process  fidelities  higher 
than  86%  are  achieved.  This  improvement  is  analogous 
to  the  early  development  in  the  superconducting  charge 
qubits  m  nu,  where  operating  at  sweet  spot  (quantro- 
nium)  with  resonant  microwaves  13  demonstrated  the 
first  high  quality  universal  single  qubit  gate  operations 
after  initial  demonstration  of  charge  qubit  manipulation 
with  non-adiabatic  pulse  techniques  [4] .  Applying  Hahn- 
echo  decoupling  sequence  provides  modest  improvement 
in  coherence  time  (T^echo  ~  2.1  ns  compared  to  T2*  ~  1.3 
ns),  indicating  that  understanding  high  frequency  charge 
noise  as  well  as  charge  relaxation  at  the  sweet  spot  will 
be  important  for  further  development.  The  quantum  dot 
charge  qubit  is  highly  tunable  using  gate  voltages,  and 
we  expect  that  investigating  coherence  and  process  fi¬ 
delity  as  a  function  of  tunnel  coupling  strength  between 
the  dots  will  provide  an  effective  route  to  improve  its 
performance  [5  . 

Methods 

Measurement:  The  experiments  are  performed  on  a 
double  quantum  dot  fabricated  in  a  Si/SiGe  heterostruc- 
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ture  E  EDI  at  base  temperature  (electron  temperature 

140  mK  [30])  in  a  dilution  refrigerator.  The  valence 
electron  occupation  of  the  double  dot  is  (2,1)  or  (1,2), 
as  confirmed  by  magnetospectroscopy  measurements  (301  • 
All  manipulation  sequences  including  microwave  bursts 
are  generated  by  a  Tektronix  70002A  arbitrary  waveform 
generator  and  are  added  to  the  dot-defining  dc  voltage 
through  a  bias  tee  (Picosecond  Pulselabs  5546-107)  be¬ 
fore  being  applied  to  gate  GR.  Similarly  to  our  previous 
study  [12],  the  conductance  change  through  the  quan¬ 
tum  point  contact  (QPC)  with  and  without  the  applied 
microwave  burst,  and  measured  with  a  lock-in  ampli¬ 
fier  (EG&G  model  7265),  is  used  to  determine  the  av¬ 
erage  charge  occupation  and  is  converted  to  the  reported 
probabilities.  Charge  relaxation  during  the  measurement 
phase  is  taken  into  account  using  the  measured  charge  re¬ 
laxation  time  T\  ~  23.5  ns  at  the  read-out  detuning  of 
—160  /ieV  [T2],  In  the  Supplementary  Information 
SI,  we  present  the  details  of  measurement  technique  and 
the  probability  normalization. 

Gate  set  tomography  (GST):  GST  is  performed  in  two 
stages,  as  described  in  Ho].  First,  data  obtained  from 
short  gate  sequences  are  analyzed  using  linear  inversion, 
to  obtain  a  rough  estimate.  Next,  data  from  long  se¬ 
quences  are  incorporated  using  maximum-likelihood  pa¬ 
rameter  estimation,  to  refine  the  preliminary  estimate. 
The  required  gate  sequences  are  defined  using  a  set  of 
fiducial  sequences  T  =  that  fix  a  consistent  (but  a 
priori  unknown)  reference  frame.  For  this  experiment, 
we  chose  T  =  {0,  Gi,  G2,  G\ },  where  0  is  the  null  opera¬ 
tion  (do  nothing  for  no  time).  In  terms  of  these  fiducials, 
linear  inversion  GST  demands  data  from  all  sequences  of 
the  form  FiGkFj  where  i,  j  =  1 ...  4  and  k  =  0 ...  3  and 
Go  =  0  is  the  null  operation.  Supplementary  Information 
S3  provides  more  details  about  the  GST  implemented 
here. 
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SUPPLEMENTARY  INFORMATION 


SI.  DETAILS  OF  TIME- AVERAGED 
MEASUREMENT  AND  PROBABILITY 
NORMALIZATION 


For  the  measurement  of  changes  in  the  probabilities  of 
charge  occupation  resulting  from  fast  microwave  bursts, 
we  use  the  general  approach  described  in  [12] ,  where  we 
measure  the  difference  between  the  QPC  conductance 


with  and  without  the  manipulation  pulse  train.  Fig.  SI 


shows  similar  scheme  adopted  in  this  work.  We  alter¬ 
nate  an  appropriate  number  of  manipulation  and  mea¬ 
surement  sequences  (measurement  time  ^60  ns)  with  mi¬ 
crowave  bursts  with  adiabatic  in  and  out  sequences  with¬ 
out  microwave  to  form  a  low  frequency  signal  with  fre¬ 
quency  on  the  order  of  ^1  kHz.  The  adiabatic  ramp  se¬ 
quence  without  microwave  manipulation  does  not  induce 
state  population  change,  but  greatly  reduces  background 
signal  due  to  capacitive  crosstalk.  The  manipulation  se¬ 
quence  including  the  microwave  bursts  is  generated  us¬ 
ing  Tektronix  AWG70002A  arbitrary  waveform  generator 
with  maximum  sample  rate  of  25  Gs/s  and  analog  band¬ 
width  of  about  13  GHz.  The  data  are  acquired  using 
a  lock-in  amplifier  with  a  reference  signal  correspond¬ 
ing  to  the  presence  and  absence  of  the  pulses,  as  shown 
schematically  by  the  orange  dashed  line  in  Fig.  [Si]  We 
compare  the  measured  signal  level  with  the  correspond¬ 
ing  (2,1)- (1,2)  charge  transition  signal  level,  calibrated 
by  sweeping  gate  GL  and  applying  the  orange  square 
pulse  shown  in  Fig.  [Si]  to  gate  GL.  Similarly  to  previ¬ 
ous  work  [12] ,  charge  relaxation  during  the  measurement 
phase  is  taken  into  account  using  the  measured  charge 


Figure  SI.  Pulse  sequence  used  for  lock-in  measurement 
of  the  qubit  state.  The  orange  dashed  line  shows  the  cor¬ 
responding  lock-in  reference  signal,  which  also  serves  as  the 
probability  calibration  pulse. 
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relaxation  time  T\  «  23.5  ns  at  the  read-out  detuning  of 
er  ~  —160  /ieV. 


S2.  MODEL  HAMILTONIAN 
CORRESPONDING  TO  FIG.  1C 


The  energy  level  diagram  is  calculated  from  the  Hamil¬ 
tonian  written  in  the  basis  of  the  ground  and  first  excited 
states  of  the  electron  charge  states  (2, 1)  =  | L)  and 

(1,2)  =  | R):  {|L)g,  |L)e,  |R)g,  \R)e}-  In  this  basis,  H  is 


given  by 


(e/2 

0 

0 

e/2  +  SEh 

Ai 

—A3 

\  —  A2 

a4 

0  -e/2  +  SEJ 


Here,  Ai_4  are  tunnel  coupling  matrix  elements,  6  is  the 
detuning,  and  SE^  and  SE  are  the  energy  separation  of 
the  ground  and  excited  | L)  and  | R)  states,  respectively. 
The  parameters  A3,  A4,  and  SE l  are  relevant  to  high  en¬ 
ergy  states  which  we  do  not  experimentally  access  in  this 
work,  and  we  use  the  values  determined  from  our  previ¬ 
ous  study  m •  We  used  parameters  of  tunnel  coupling 
between  | L)  and  \R)  2Ai  jh  =  5.2  GHz,  tunnel  coupling 
between  |  L)  and  low  lying  excited  state  |  R)e  on  the  right 
dot  2A2 /h  =  14.5  GHz,  and  singlet-triplet  energy  split¬ 
ting  on  the  right  dot  SR  =  12.1  GHz,  extracted  from 
measurements  as  in  Refs.  HU  HU. 


S3.  STATE  AND  PROCESS  TOMOGRAPHY 

A.  Details  of  the  linear-inversion  gate  set 
tomography  experiments 

Gate  set  tomography  (GST)  is  a  method  for  self- 
consistently  characterizing  state  preparation,  measure¬ 
ment,  and  quantum  operations.  Our  method  is  described 
in  detail  in  Ref.  [TO]  (Ref.  [3Tj  demonstrated  a  similar 
protocol);  this  section  discusses  some  specific  issues  en¬ 
countered  in  this  experiment. 

Like  standard  QPT,  GST  requires  an  informationally 
complete  set  of  measurements  -  i.e.,  at  least  one  linearly 
independent  probability  per  gateset  parameter.  But  it 
eliminates  standard  QPT’s  need  for  known  initialization 
and  measurement  processes  by  doing  initialization  and 
measurement  with  unknown  operations  selected  from  the 
gateset  to  be  characterized.  Informational  completeness 
is  achieved  by  using  fiducial  sequences  (F  =  {F^})  built 
from  the  gates  themselves  to  generate  at  least  d2  =  4 
distinct  input  states,  and  to  rotate  the  final  state  into 
at  least  d2  =  4  distinct  measurement  bases.  However, 
whereas  in  standard  QPT  informational  completeness 
can  be  ensured  in  advance,  in  GST  we  do  not  know  in 
advance  what  operations  the  Fi  will  perform.  Instead, 
we  test  the  data  for  the  signature  of  informational  com¬ 
pleteness,  and  if  necessary  adjust  the  fiducial  sequences. 


This  signature  is  the  spectrum  of  the  Gram  matrix  G 
defined  by  Gij  =  Pi(FiFj).  G  should  have  rank  d2 
if  the  fiducials  T  are  informationally  complete  for  a  d- 
dimensional  system.  In  the  experiments  reporte  here,  our 
initial  choice  of  T  accidentally  produced  Gram  matrix 
with  only  three  significant  eigenvalues;  we  adjusted  the 
fiducial  sequences  to  ensure  that  all  four  eigenvalues  of 
G  were  large  (>  0.19).  With  our  initial  microwave-pulse 
gates,  chosen  to  roughly  approximate  7r/2  rotations,  the 
fiducial  set  T  =  {0,  Gi,  G2,  Gi2}  yielded  a  Gram  matrix 
with  smallest  eigenvalue  of  0.05,  which  caused  the  linear 
inversion  GST  to  be  ill-conditioned.  Changing  the  fourth 
fiducial  to  F4  =  Gf  as  well  as  reducing  the  rotation  an¬ 
gle  by  changing  the  amplitude  of  the  ac  driving  yielded 
a  Gram  matrix  with  smallest  eigenvalue  0.19,  which  en¬ 
sures  reliable  linear  inversion  for  the  first  step  of  GST. 
As  is  shown  in  Fig.  4c  of  the  main  text,  GST  estimates 
the  rotation  angles  of  Gi,  G2,  and  G3  to  be  0.387T,  0.417T, 
and  0.457T,  respectively. 


B.  Details  of  the  non-linear  gate  set  tomography 
refinement 

GST  proceeds  in  two  steps:  a  linear  inversion  step, 
where  we  obtain  coarse  estimates  of  the  gates,  followed 
by  a  refinement  step.  Higher  accuracy  is  attained  by  per¬ 
forming  longer  sequences  of  gates.  In  this  experiment, 
we  performed  (and  estimated  the  measurement  probabil¬ 
ity  associated  with)  sequences  of  the  form  F^GJiFj,  for 
i,  j  m  1 . . .  4,  k  =  1 . . .  3,  and  n  =  2  . . .  32.  Many  repeti¬ 
tions  amplify  over/under-rotation  error  (e.g.  by  an  angle 
6)  in  each  gate  G/c,  so  that  by  measuring  nO  to  reasonable 
accuracy  we  achieve  a  very  accurate  estimate  of  9.  How¬ 
ever,  because  the  associated  measurement  probabilities 
are  highly  nonlinear  in  G&,  we  cannot  use  linear  inver¬ 
sion.  Instead,  we  estimated  the  gateset  Q  =  {G/c,p, M} 
by  maximizing  the  likelihood  C{Q)  =  Pr(data|t/).  Since 
the  likelihood  function  C{Q)  is  in  general  nonconvex,  we 
used  the  linear  inversion  estimate  as  a  starting  point 
for  local  optimization.  We  found  that  in  this  case  the 
best  local  optimization  procedure  (resulting  in  the  high¬ 
est  maximum  likelihood)  was  to  repeatedly  apply  the 
Nelder-Mead  downhill  simplex  method  until  it  converged. 
Since  the  resulting  estimate  is  not  necessarily  a  physically 
valid  gateset  (those  for  which  p  and  M.  are  valid  quan¬ 
tum  states/measurements,  and  each  G&  is  a  completely 
positive,  trace-preserving  map),  it  is  projected  onto  the 
space  of  valid  gatesets  to  produce  a  final  estimate.  This 
final  maximum-likelihood  estimate  predicts  count  statis¬ 
tics  that  fit  the  observed  data  well,  but  the  variability  of 
results  with  different  local  optimization  procedures  indi¬ 
cates  that  the  likelihood  contains  a  number  of  local  min¬ 
ima  and  suggests  that  even  in  our  best  case  the  likelihood 
may  not  have  been  strictly  maximized.  However,  this  ef¬ 
fect  appears  to  be  less  significant  than  non-Markovian 
noise  effects,  which  are  not  accounted  for  in  GST  (unlike 
QPT,  GST  can  incorporate  data  from  experimental  se- 


quences  in  which  a  single  gate  Gk  appears  more  than  once 
-  but  it  does  so  by  explicitly  assuming  that  Gk  applies 
the  same  quantum  process  each  time  it  is  applied,  with 
no  systematic  variation  with  respect  to  time  or  context). 


C.  Standard  quantum  process  tomography  (QPT) 

We  perform  standard  process  tomography  on  our  qubit 
in  the  usual  way,  by  preparing  precalibrated,  informa¬ 
tionally  complete  states,  applying  the  unknown  operation 
to  be  characterized,  and  measuring  in  a  precalibrated, 
informationally  complete  basis.  In  our  experiment,  the 
state  preparation  and  measurement  processes  are  precal¬ 
ibrated  using  Rabi  oscillation  and  Ramsey  fringe  experi¬ 
ments,  as  detailed  in  the  main  text.  We  perform  QPT  to 
characterize  a  total  of  four  processes:  X^/2,  Z^/2  (shown 
in  Fig.  4a  in  the  main  text),  and  G\  and  G3  (shown  in 
Fig.  4b  in  the  main  text).  In  all  cases,  first  we  initial¬ 
ize  to  the  | L)  state,  then  adiabatically  transition  to  the 
|0)  state,  as  detailed  in  the  main  text.  Next,  we  prepare 
an  informationally  complete  set  of  inputs  in  two  sepa¬ 
rate  ways.  First,  for  the  X^/2  and  G\  (which  is  X-like) 
operations,  we  apply  the  operations  X^/2Z^/2  and  X$ 
for  three  separate  angles  4>.  This  gives  the  \—X)  state, 
in  addition  to  three  states  that  he  on  the  y  —  z  plane, 
forming  an  informationally  complete  set.  Second,  for  the 
Z^j2  and  Gs  (which  is  Z-like)  operations,  we  either  do 
nothing  or  apply  the  operations  Z^Xn/ 2  for  three  sep¬ 


arate  angles  <F.  This  gives  for  our  inputs  the  states  |0) 
and  three  states  in  the  x  —  y  plane,  which  are  again  in¬ 
formationally  complete.  The  reference  time  which  sets  <f> 
is  varied  in  the  experiment  and  the  results  with  10  differ¬ 
ent  <f>  were  used  to  estimate  statistical  average  and  stan¬ 
dard  deviation  of  the  process  matrices  and  fideltities.  We 
choose  these  two  different  input  sets  for  X  or  Z-gate  pro¬ 
cess  tomography  for  data  post-processing  convenience. 
After  preparing  the  informationally  complete  input  set, 
we  then  implement  the  processes  to  be  characterized.  Fi¬ 
nally,  we  perform  —Yn/ 2,  X^/2,  and  identity  operations, 
followed  by  a  measurement  in  the  Z-basis.  Note  that  in 
the  state  tomography  of  X-gate  we  implement  tomogra¬ 
phy  in  the  rotating  frame  by  adjusting  phase  of  the  mea¬ 
surement  7r/2  pulses  taking  phase  accumulation  during 
manipulation  pulse  into  account,  whereas  Z-gate  tomog¬ 
raphy  is  performed  in  the  lab  frame  where  the  phase  of 
the  measurement  pulses  are  fixed.  Once  these  measure¬ 
ments  are  complete,  we  use  maximum-likelihood  estima¬ 
tion  00  El  to  reconstruct  the  x  matrix  representation 
of  the  unknown  process,  given  by  m 

4 

£{p)=  X  ZmpElxmn ,  (S2) 

m,n=l 

where  £(p)  is  the  density  matrix  specifying  the  output 
for  a  given  input  density  matrix  p,  Em  are  the  basis  op¬ 
erators  in  the  space  of  2  x  2  matrices,  and  x  is  the  process 
matrix.  The  entries  of  this  matrix  are  plotted  in  Fig. 4a 
of  the  main  text. 
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